Dilepton production from hot hadronic matter in nonequilibrium 



B. SchenkcQ and C. GreineiQ 

Institut fur Theoretische Physik, Johann Wolfgang Goethe - Universitat Frankfurt, 
Max-von-Laue-Strafie 1, D-60438 Frankfurt am Main, Germany 

The influence of time dependent medium modifications of low mass vector mesons on dilepton 
yields is investigated within a nonequilibrium quantum field theoretical description on the basis of 
the Kadanoff-Baym equations. Time scales for the adaption of the spectral properties to changing 
self energies are given and, under use of a model for the fireball evolution, nonequilibrium dilepton 
yields from the decay of p- and a;-mesons are calculated. In a comparison of these yields with 
those from calculations that assume instantaneous (Markovian) adaption to the changing medium, 
quantum mechanical memory effects turn out to be important. 



PACS numbers: 11.10. Wx;05.70.Ln;25.75.-q 

Keywords: nonequilibrium quantum field theory; relativistic heavy-ion collisions; dilepton production 



I. INTRODUCTION AND MOTIVATION 

Relativistic heavy ion reactions, as performed at the SIS at GSI, Darmstadt, the AGS at BNL, the SpS at CERN 
and the RHIC at BNL, allow for studying strongly interacting matter under extreme conditions at high densities and 
temperatures. One of the main objectives is the creation and identification of new states of matter, most notably 
the quark-gluon plasma (QGP). Photons and dileptons do not undergo strong final state interactions and thus carry 
undistorted information especially on the early hot and dense phases of the fireball. Photon spectra are a suitable 
observable for the temperature of such a system whereas dileptons have encoded additional dynamic information 
via their invariant mass. Particularly in the low mass region, dileptons couple directly to the light vector mesons 
and reflect their mass distribution. They are thus considered the prime observable in studying mass (de-) generation 
related to restoration of the spontaneously broken chiral symmetry. Additionally,vector mesons are also affected by 
many-body effects due to coupling to baryonic resonances. For an overview see m. Indeed the CERES experiment 
at the SPS at CERN 0, has found a significant enhancement of lepton pairs for invariant masses below the pole 
mass of the p- meson, giving evidence for such modifications. 

In order to be able to extract precise information from the data, it is essential to find a thorough mathematical 
description for the dilepton production of an evolving fireball of strongly interacting matter. In particular it is 
necessary to consider the fact that especially during the early stages of a heavy ion reaction the system is out of 
equilibrium - the medium and hence the properties of the regarded mesons undergo substantial changes over time. 
These scenarios have been described within Boltzmann-type transport calculations using some quantum mechanically 
inspired off-shell propagation in [3, IH Itl l?l IsL Hcllllj . I n principle, a consistent formulation beyond the standard 
quasi-particle approximation is needed - fully ab initio investigations of off-shell mesonic dilepton production without 
any further approximation does not exist so far. 

More generally nonequilibrium quantum field theory has become a major topic of research for describing microscopic 
transport processes in various areas of physics (see e.g . |T^ ) and references therein). One major question deals with 
how quantum systems eventually thermalize. In [12j the quantum time evolution of </> 4 -theory for homogeneous 
systems in 2+1 space-time dimensions for far from equilibrium initial conditions has been investigated, while earlier 
works (e.g. 0) studied the 1+1 dimensional case. It was shown in 01 that the asymptotic state in the far future 
corresponds to the exact off-shell thermal state of the system obeying the equilibrium Kubo-Martin-Schwinger (KMS) 
relations among the various two-point functions. For a coupled fermion-boson Yukawa-type system in 3+1 dimensions, 
eventual equilibration and thermalization in momentum occupation was shown in jbl| . In addition, the full quantum 
dynamics of the spectral information was analyzed in [T^ . This issue was also addressed in some detail in |15). In 
a subsequent work |16| the exact solutions were examined in comparison to approximated, instantaneous off-shell 
transport equations, obtained by a first order gradient expansion. In that particular case it turned out that indeed 
these approximated equations are a very good substitute for the full dynamics. 
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When dealing with vector mesons the important question emerges whether such a quasi instantaneous adaption of 
the dynamic and spectral information to the changing medium, as also assumed in more schematic model calculations 
0, an d Monte-Carlo kinetic transport simulations ^3], is always a suitable assumption. In a ultrarelativistic 
heavy ion collision the typical lifetime of the diluting hadronic phase is only 5 — 8 fm/c 1]. On the other hand the 
spectral information should always react on temporal changes with a certain " quantum mechanical" retardation. If the 
timescale of these changes of the system becomes comparable to the retardation time, an instantaneous approximation 
becomes invalid and memory effects for the spectral properties of the excitations are present. Such possible (quantum 
mechanical) memory effects, i.e., potential non-Markovian dynamics, have often appeared in descriptions of the 
microscopic evolution of com plex qu antum systems. Also a nonequilibrium treatment of photon production from a 
hot QGP was given in [l9l |20L l2lL l23 | and the importance of memory effects and shortcomings of the S-matrix approach 
were pointed out for that case. The question of whether memory effects arc important for dilepton production from 
hot hadronic matter, i.e., whether such effects have influence on the measured dilepton yields, constitutes the major 
motivation for the present study. 

We give for the first time explicit calculations for dilepton production from first principle nonequilibrium transport 
equations. We approach the problem using a nonequilibrium quantum field theoretical description based on the 
formalism established by Schwinger and Keldysh |23l |2J, l25l 26]. By this we set up a framework incorporating 
the full quantum dynamics, which is necessary for the description of the transport of true off-shell excitations. We 
present a new derivation of the formula for the dynamic dilepton production rate starting from the Kadanoff-Baym 
equations |27|. which are nonlocal in time and hence account for the finite memory of the system. The resulting 
formula for the rate involves a (half) Fourier transform over past times of the two-time Green function of the virtual 
photon. The issue within this representation is that the rate for an invariant mass in the range of interest is hidden 
as a tiny component in the two-time function, which contains the rate for all invariant masses. With this work we 
meet the challenge of evaluating this expression for the nonequilibrium dilepton production rate, which is particularly 
important since it is the only causal approach that retains all memory effects. Any treatment so far, involving gradient 
expansions of the Kadanoff-Baym equations, where future contributions to the Green function are treated equally to 
those from the past, can not precisely describe a system that is quickly evolving with respect to the timescales on 
that the regarded quantities adjust to system changes. The reader, who is already familiar with the formalism of 
nonequilibrium quantum field theory may skip most of the first part of Section ITU and continue reading with Eo. (|I9|l . 

The paper is organized as follows. We start with a brief introduction of the used formalism and the presentation of 
our new derivation of the dilepton production rate for nonequilibrium systems in Sections III Al and III Bl The involved 
two-time Green function of the virtual photon is further discussed in Section III CI It is then shown how the medium 
modifications of the light vector mesons enter the dilepton rate via the principle of vector meson dominance. The 
simulation of medium modifications of these vector mesons by introduction of a certain time dependence of their self 
energies is introduced in Section IlII Al We analyze the contributions to the rate in time representation and discuss 
the precision of the numerics in Section IIII Bl A quantitative description of the retardation is given in Section IIII CI 
by introduction of time scales, which characterize the memory of the spectral function. Comparing them to typical 
time scales in heavy ion reactions reveals that changes of the spectral function can not generally be assumed to 
be adiabatic and that memory effects can become important. We discuss quantum mechanical interference effects 
occurring within this full quantum field theoretical description in Section IlII Dl Dilepton yields are calculated first 
for constant temperature and volume in Section fill El Finally, convolution of the dynamic rates with a fireball model 
employing a Bjorken like expansion leads to our most important result, presented in Section IlII Fl Comparison to 
the quantities computed in the static limit, where all meson properties adjust to the medium instantaneously, the so 
called Markov limit, reveals the significance of memory effects and the consideration of the full dynamics for certain 
cases such as the celebrated and continuously discussed Brown- Rho scaling |28| . 



II. THE NONEQUILIBRIUM PRODUCTION RATE 
A. Lepton number transport equation 

We utilize the Schwinger-Keldysh realtime formalism and the emerging Kadanoff-Baym equations [27| in order to 
derive the dynamic nonequilibrium rate of produced electron-positron pairs, coming from the decay of light vector 
mesons via virtual photons in a spatially homogeneous, yet time dependent system. A different derivation for the 
dilepton rate was performed in .29] for dileptons from a pion plasma as well as in |30j. starting from the dilepton 
correlator. The resulting formulas will provide a powerful tool to compute the dynamic behavior of the dilepton 
production rate, influenced by a changing surrounding medium. 

We extract the number of produced electrons with momentum p at time r from the Wigner transform of the 
electron propagator G < (1, 2) = «(4 , (2)5'(1)) at equal times t\ = t<i = r, which for a general system is given by (using 
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notation as in pjll]') 

v P+ s> P+ ., s ) «(p + 1, S )«(-p + 1, r)e ^x e ir ( - E+ - B _) j i (1) 

by projecting on the quantity (& pr 6 P)S ) <5 rs |3^ |: 

N(p,r)=-if d 3 XTr {V G< (X, p, r) } . 
It can be easily verified that this is achieved by use of the projector 

V = Jojr^2u(P,sW (p,s) = 7o — L(|> + to)7 , (2) 

where we used J^s u (Pi s)u(p, s) = + to). £+ = -Ep+a and = -E p _a are the energies corresponding to 

momentum states ±(p + § ) and ±(p — ?f). For further details on the spin decomposition of the Wigner function see 
jHHHH The equations of motion for G<(1,2) and G£Yl,2) = -i(*(l)*(2)) are the Kadanoff-Baym equations, 
generalized to the relativistic Dirac structure 32, 33, 34, 35]: 



(iTX-m-Sjj^ljJG? (!,!')= / d2(E>(l,2) - E<(1, 2))G*(2 > 1') 



d2S^(l,2)(G>(2,lO-G<(2,l')), (3) 



1') (-fy^, - m - Shf(1')) = j[ 1 d2(G>(l, 2) - G<(1, 2))E*(2, l') 

- r'd2G^(l,2)(E>(2,l')-E<(2,l')), (4) 

•/to 

with the self energy E and its local, Hartree-like term Hjjf- (1,2) is the short term notation for the coordinates 
(t\, Xi, i 2 , x 2 ). Using retarded and advanced Green functions 

G ret / adv (l,2) = ±6(±(h -t 2 ))(G>(l,2)-G < (l,2)), (5) 
an important relation can be obtained directly from the Kadanoff-Baym equations: 

G^(l,l')=/ d2 / d3G rot (l,2)E^(2,3)G adv (3,l') 

+ j dx. 2 j dx 3 G rct (l,x 2 ,i )G^(x 2 ,io,x3,to)G adv (x 3 ,to,l')- (6) 

It can be regarded as a generalized fluctuation dissipation relation j3j 13^. l3^ | . The second term accounts for the 
initial conditions at time to only. It can be neglected if one lets the system evolve into a specified initial state for a 
sufficiently long time. This is done by keeping the self energy insertions time independent prior to the onset of the 
dynamics at the initial time t$. Fourier transformation of @ and J3J in the spatial coordinates (xj — xi<) and taking 
7o© — 07o at t\ = ty = t leads to 

id T G< (p, r) - p • ( 7 o7G< (p, r) - G< (p, t) 77o ) 

- m( 7o G<(p, r) - G< (p, r) 7o ) = 7o G (p, r) - G (p, r) 7o , (7) 
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with the collision-terms 

G(p,r) = f df(E>(p,r,f)G < (p,f,T)-E<(p,r,i)G > (p,f,r)) 
J to 

C(p,r) = f di(G > (p,T,t)E<(p,i,r) - G < (p,T,t)E>(p,i,T)) , 

J to 

where Hhf has been effectively absorbed into the mass m. Application of the projector © to Eq. Q then yields 
the electron production rate at time r: 



d T N(p,r) = -Tr{p( 7o G(p,r)- G(p,r) 7o )} 
= (-2)Rc[Tr{7'(7oG(p,T))}" , 



(8) 



Due to having a very long mean free path, the electrons are not expected to interact with the medium after they have 
been produced. This is why they can be described using the free propagators 



G£(p,t,T) = iT^rhoEp + j ■ p - m)e 



iEp(t-r) 



2R 



G > (p,i,r) 



2E X 



■(loEp -7-p + m)e 



iE p (t-T) 



(9) 
(10) 



With that Eq. JSJ) becomes 

d T N(p, t) = 2Re 
= 2Rc 



where in the second step we used that Gq (p, t, t)V^q — Gq (p, t,r) and G^(p, t, t)V^iq — together with the cyclic 
invariance of the trace. Inserting (|10|> finally leads to 



Tr jp 7 o Qf dt (E< (p, r, i)G > (p, i, r) - E> (p, r, t)G < (p, i, r)) 
Trj^ ^(S < (p,r,i)G > (p,i») 





r r^ + m r 


2 Im 






L l 2 ^p Jto 



(11) 



with po = Ep- The full dynamic information for the production of an electron at a given time r is incorporated in 
the memory integral from the initial time to until the present r on the right hand side of Eq. ■ 



B. The electron self energy E 

The medium as the source for the production of dilcptons enters via the dressing of the virtual photon propagator 
in the electron self energy (see Fig. This dressing will finally be given by vector mesons. 




FIG. 1: Feynman graph for the electron self energy E(l, 2) 
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II is the self energy of the virtual photon and we have 



^I] < (p,^l,^2) - :J - 



(12) 



with _D5" A " y (k, t\ 1 t-i) being the propagator of the dressed virtual photon with momentum k. Inserting this with the 
explicit form of the free electron propagator @ into Eq. (jl 1|> yields 



d T N(p,r) = 2 Re 



dt 



to 



d 3 k 

(2^)3 



(i)D<^(k, r,f)e iE p( T -*)e iS p- k ( T -*) 



x — — Tr 

2E p 2_E p _k 



i)7/*(7D- E P-k + 7 • (P - k ) - 



(13) 



where evaluation of the trace leads to 



d 3 k 1 1 



) 3 E p £ p _ k 



[Pp.{k - p) v +p v {k- p)^ - g^u(Pv,(k - pY + m 2 )] 



x Re 



f dtiD<^(k,T,i)e l(E *> +Ek -*> ) 

Jtn 



(14) 



i(E p +E k - p )(T-t) 

Defining p~ = k — p and p + = p as the four-momenta of the outgoing electron and positron, we rewrite Eq. (|14|) to 

E+E. 



( r ) =7^3b [P/fo- + PZPn ~ 9^(p + P + ' 



d 3 p + d 3 p 



(2tt) 6 l ^ 
x Re 



f diiD<^(k,T,t)e 

J to 



i(E++E-)(r-t) 



(15) 



with _E_ = i?k-p and E + = E p . R denotes the number of lepton pairs per unit four- volume, produced with the 
specified momentum configuration. We now show that Eq. (|15f) is indeed the generalization of the well known 
thermal production rate for lepton pairs [38ll39ll4Cl| . Using the Fourier transform in relative time coordinates, defined 
by 



D<^(k,r,t) = J ^-D<^(k,r,u;) t 



— itu(r— i) 



(16) 



and taking t — > — oo, we have for the stationary case: 



dR 



2e 2 



E + E - i • ,.>_',,; L J '/' 



d 3 p+d 3 p~ (2tt) 



-9»v(p + P +m 2 )] 



x Re 



7 27r 7_„ 



(17) 



^< M !/ (l C) (jj) i s rea j an( j time independent in the stationary case. The real part of the last integral is simply ttS(E^ 
E- — ui). With that and the virtual photon momentum fc M = (E, k), £7 = E + + E-, the rate becomes 



dR 



+ d 3 p+d 3 p 



(27T) 



[P>; - + m 2 )] £><^(fc) 



1 



(2tt) 6 j 
2e 2 



[ptvZ +PZPZ. - 9AP + P~ + ™ 2 )] T74 n 7 ""W 



(2tt) 6 



P^Pis +PZP» ~ 9nv(p + P +m 2 )] 



M 4 7 
1 1 



imir ct ^(fc), 



(18) 



where we used Eq. © for D< ^(k) in its stationary limit in the first step as well as n < = 2mBlmII ret , which follows 
from the Kubo-Martin-Schwinger (KMS) relation (37L EH |42 | . in the second step. Eq. i|18|) is the well known rate of 
dilepton production derived in e.g. |40| 1 . 



1 The difference in the overall sign is due to the opposite sign in the definition of the Green functions in |4(1 . 
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We return to the nonequilibrium formula (|15|l and project on the virtual photon momentum using 



dN 



dR 



6 4 {p + +p~ -fc)dV(fp 



This leads to 



-(r, k, E) 



d 4 xd 4 k J d 3 p+d 3 p 
d 3 p + f d 3 p 



d 4 xd 4 k' 



(2tt) 6 
x Re 



P_ 



6 4 {p + +p - k)[p+p u +ptp^ - 9^(P + P +m 2 )] 



f dtiD<^(k, T ,t)e 



iE(r-i) 



(19) 



for the production rate of dilepton pairs of momentum k = (E = E + + P_ , k) . In the following numerical study we 
will consider the mode k = exclusively, i.e., the virtual photon resting with respect to the medium. For this case, 
and taking the electron mass m to zero, Eq. 119|) can be simplified to 



dN . , „ rt 2e 2 2 , , 

-ir.k 0,E) =- ' ! - '• '- 



d 4 xd 4 k 



(2tt) 6 3 



^{k^ky - k g^Re 



f dtiD< (k = 0, t, t)e iB ( T -^ 



(20) 



This expression is easily understood. The dynamic information is inherent in the memory integral on the right hand 
side that runs over all virtual photon occupation numbers, Fourier transformed at energy E from the initial time to 
the present. Hence this memory integral determines the full nonequilibrium dilepton production rate at time r. 



C. The in-medium virtual photon self energy n 

The dilepton production rate 1121 III involves the virtual photon occupation number, expressed by the propagator 
D< . We introduce the dynamic medium dependence by dressing this virtual photon propagator with the medium 
dependent p- or cj-meson. This dressing enters via the photon self energy II < through the fluctuation dissipation 
relation (cf. ®) for D<» v : 

D< = Pf 1 a n<3 Pf (21) 

where implies the integration over intermediate space-time coordinates. In the medium the vector mesons and 
virtual photons have two possible polarizations relative to their momentum in the medium. This leads to two different 
self energies IIt (transverse) and IIl (longitudinal). Introduction of the projectors Pl and Pr allows us to split the 
propagators and the self energy into a 3-longitudinal and a 3-transversal part, relatively to the particle's momentum 

pfia pfj,a p/3u p/3v 
jjTetfia _ ■'T L . jj&dv (3v _ r T r L zoo") 

7 ~ k 2 - n** k 2 - njf ' 7 k 2 - n!f* k 2 - n™** 1 ' 



and 



Kb = ~Pt a pB% - P L a/3 U< . (23) 



The projectors fulfill the usual projector properties P(t/l) = ^(t/l) an d PtPl = PlPt = 0. With that, Eq. I|21|) 
becomes 



d<»» = pf ^ (-p T a/3n< - P w n<) Df v ^ 

7,T T ^j.L^h ' 

with 



= -D<P^-D<P^, (24) 



^T/L ~ k 2_ nrc^ © n T/L © fc2 _ jjretj. ■ (25) 

In the case k = 0, considered later, the longitudinal and transverse parts become identical and it follows 



n <^__ n < piiv _ n < p^ _ _ D < I pliv , piiv\ _ _n< „, 

u 1 — U 7,T r T 7.T L — ^,Tr T ' L / 7,T \9 ^ 2 
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In this case the production rate depends only on the transverse part of the virtual photon propagator. The appearing 
factor of 3 accounts for the two transverse and one longitudinal directions: 



/ diiD< T {\L = Q,T,?)e iB{ ' T - h 

J t 



(26) 



For dilepton production II rct oc e 2 and E is the invariant mass of the virtual photon. For the cases, we are interested 
in, it holds |II ret | <C E and we can approximate 

D< T = D; rt o 0n<0D*. (27) 

For k -> 

D$(k -> 0, t) = -0(t)t = D^(k -> 0, -t) (28) 
and we may calculate the virtual photon propagator using the transport equation 

»D< T (k = 0,T,t)= / dti / dt 2 (r-i 1 )(*II$(k = 0,ii,t2))(t-t2). (29) 

•/to "/to 

The appearing undamped photon propagators lead to diverging contributions from early times, i.e., for low frequencies. 
In the later numerical calculation these contributions turn out to be at least of the order 10 5 larger than the actual 
(higher frequency) structure. Due to the naturally limited numerical accuracy, the higher frequency structure would 
get lost among these early time contributions. In order to cure this numerical problem, we introduce an additional 
cutoff A for the free photon propagators, i.e., we perform the replacement: 

Dy%r - h) = (r- h) -f (r - ti)e- A < T -*0 (30) 

and analogously for D^fa — t). In the performed calculations we employ A as 0.3 GeV. The exponential factors lead 

to a reduction of the rate, which we will overcome by renormalizing the final result by multiplication with — ^ — 
a factor we get from Fourier transforming the convolution (|29|) . assuming equilibrium. This will not affect the time 
scales we are interested in, and comparison of the dynamically computed rate for the stationary case (constant self 
energy) with the analytic thermal rate shows perfect agreement. 

Vector meson dominance (VMD) HI Q allows for the calculation of 11^ , using the identity between the electro- 
magnetic current and the canonical interpolating fields of the vector mesons |45|: 



e 



Jn = m P Pn (31) 

9p 



which leads to 



2 



Bp 



for the self energy. When treating the w-meson, we use the corresponding self energy and propagator. We again apply 
the generalized fluctuation dissipation relation 10 to calculate 

fl< T = ^0S< T 0fl* (33) 

with the p-meson self energy T . The transverse parts of the retarded and advanced propagators Dp C ^(k, t\, £2) = 
_Dp d r^(k, t2, ii) of the vector meson in a spatially homogeneous and isotropic medium follow the equation of motion 

(~0l -ml- k 2 ) D^(k, h,t 2 ) - t dtS^(k, h,t)Df T (k, t, t 2 ) = *(t x - t 2 ). (34) 

Jt 2 

In the following we will omit the index T for convenience. 

The dynamic medium evolution is now introduced by hand via a specified time dependent retarded meson self 
energy S rot (r, u>) with system time r (see Section fill Afl . From that the self energy S < , needed for solving Eq. I|33|l . 
follows by introduction of an assumed background temperature of the fireball. The fireball, constituting the medium, 
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generates the time dependent self energy E rct and, assuming a quasi thermalized system, the p-meson current-current 
correlator E < is given via 

E < (r,w,k) = 2m B (r(r))ImE rct (r,w,k), (35) 

which follows from the KMS relation 

E<(w,k) = Te" /3 "S > (^,k), (36) 

being valid for thermal systems [^3, El E3 ■ n B is the Bose-distribution. The assumption of a quasi thermalized 
background medium is of course rather strong, but necessary in order to proceed: In principle, for a full nonequilibrium 
situation the self en ergi es E < and E> for the p-meson have to be obtained self-consistently via e.g. coupling to 
resonance-hole pairs [46|, being out of equilibrium themselves. (For a realization of true nonequilibrium dynamics of a 
homogeneous system within a <fr 4 -theory see |T^ and within a coupled fermion-meson system |l4j.) For an expanding 
and inhomogeneous reaction geometry this is still not possible today. In any case, explicit calculation of E < in the 
two-time representation will cause even stronger memory effects. Additionally, in order to simulate a realistic situation 
by application of a fireball model, a temperature needs to be defined and hence local equilibrium has to be assumed. 

The framework, we have established, allows us to calculate the dynamic evolution of vector mesons' spectral 
properties, occupation, and the nonequilibrium dilcpton production rate, given an evolving medium with assumed 
time dependent quasi-thermal properties. At this point it is also worthwhile to point out that working in the two-time 
representation has great advantages over the mixed representation, in which all quantities are expressed by their 
Wigncr transforms. The problem in this case is that the Wigner transforms of the many appearing convolutions are 
nontrivial. They can be expressed by a gradient expansion to infinite order |47^ . Explicit calculations are usually 
carried out within a first order gradient approximation, which is not applicable in our case because the system is not 
evolving slowly with respect to the relevant time scales (see Sections IIII CI and IIII Fl for details on the time scales 
involved). All memory is lost by application of this approximation, because when using the mixed representation, full 
Fourier transformations from relative time to its conjugate frequency are involved, which imply treating contributions 
from the future and from the past on equal grounds. This violates causality in a quickly evolving system. Hence, 
calculations within the two-time representation allow for the most exact investigation of the evolving system under 
consideration. 



III. NONEQUILIBRIUM DILEPTON PRODUCTION FROM AN EVOLVING MEDIUM 



A. Vector meson self energies 



We are now able to calculate the evolution of the spectral properties of vector mesons in a changing medium as 
well as the corresponding dilepton production rates. We will treat mass shifts described by Brown- Rho scaling |28| . 
broadening, as caused by pion scattering, and scattering of the mesons with nucleons, leading to further broadening 
and excitation of resonances 0, 0, 0, Hfl 001 E3> E3 • The medium effects are introduced via a specific time evolving 
self energy for the vector mesons. The main purpose of this work is to investigate medium modifications dynamically 
for the first time, and compare the results to those obtained from instantaneous, Markovian calculations. 

In the following we first employ simplified self energies for possible broadening such as 

ImE rct (r,w) = -wr(r), (37) 

with a k- and w-independent width T, which, in the static case, leads to a Breit- Wigner distribution of the spectral 
function 

A(cj,k) = -hmD^(cj,k) = \ ^ (38) 

7T H 7T \uj — k — vn 2 ) 2 + (ujr) 2 

The self energy Ij37f) is given in a mixed time-frequency representation. The time dependence is being accounted for 
by introduction of the system time r, for which we will discuss two possible choices below, to stems from the Fourier 
transformation in relative time (t\ —tz). 

Modifications of the mass are introduced directly as via a local term in the self energy. From this sort of self energy 
several numerical issues emerge. In Eq. I|34[l the self energy enters in time representation 

£ rot (ii,i 2 ) = t 2 ) - ^ < (ti,t 2 ))9(t 1 -t 2 ), (39) 
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Its Fourier transform in relative time (t\ — 12), E ret (r, ui), is given by a convolution of the ^-function's Fourier transform 
and (E> - E<)(t,w): 



S ret (r,w) = i 



1 



VP 



2ir uj — ui - 
du> 1 



: (E>-E<)(t,u;) 



where we used the standard relation 



2ir ui 



1 



With (E> - S < )(r,w) 



cj — w + «e 
-2iwr(r), one has 



E<)(r,a)) + -(E>-E<)(r jW ) 



— iir5(u> — cj). 



(40) 



(41) 



E rct (r,^) = 27> f^— r(r)- !W r(T), 

' 27T CJ — UI 



which has the imaginary part (|37(l . In addition we also get an "unwanted" dispersive real part that causes an (infinite) 
mass shift. The introduced cutoff c cures that infinity and we can renormalize the mass by the replacement 



y/m 2 - ReE ret (w,r), 



with 



ReE ret (r,w) = Re 



2T(t)V 



+c 



dui u> 



2tt ui 



2 . , ( ( 1 I uj 2 



,4 \ 



(42) 



with the hypergeometric function 2-Fi, Taylor expanded to second order in the last line. To make the renormalization 
^-independent, we replace ui by the physical pole mass m, causing only minimal inaccuracies far from the peak 
position. 

Another potential problem is the behavior of the self energy E < (r, ui) for negative frequencies. As already described, 
for a thermalized system, it is given by E< (r, ui) = 2in^{T(r) 1 w)ImE Ict (T, to). For negative frequencies the propagator 
D K (r, us) contains a factor of (1 + n(ui)) since there holds the symmetry relation D K (r, — io) — D > (r, oj) for the scalar 
boson case and D > (t, u>) = 2i(l + tib)A(t, ui) in equilibrium. The first term in parentheses leads to divergent vacuum 
contributions, which are uninteresting for the investigation of the positive frequency behavior and cause numerical 
problems. To see this, we first give the limits of the expression E < (r, ui) in frequency representation: 



oc e -"/ T ( r ) for uj > 
iE < (r,w) = { 2T{t)Y{t) for ui = 
-2wr(r) for u < 



(43) 



The Fourier transform thus becomes a ^'-distribution for the negative frequency part. In order to cure this problem 
we first split the self energy into positive and negative frequency parts 



iE<(r,w) 



e M/T(r) _ I 
= :iS<. rff (r,w) 



-2uiT{t)6{-w), 



(44) 



and neglect the negative frequency part, violating the mentioned symmetry relation for D < and D > (also for E < and 
E > ). This violation however only means an omittance of the vacuum contributions for negative frequencies and does 
not cause changes for positive frequencies. The same approach may be taken with other types of self energies having 
the same symmetry property under the exchange u> «-> —ui. 
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energy for coupling to J p = | -resonances is given by jS^I (also see [54| , who treated coupling of pions to resonance 



Coupling of resonances to the p — TV-channel has been treated by 0, lEH l53| . For the k — mode, the full self 
lergy for coupli] 
hole excitations) 

lu 3 ET r (t) 



ImS(r ) w,k= 0) 



p(t) ( fRNp \ 

3 V m p J 1 



(w 2 - E 2 



(45) 



E = 



i% + k — in 7v is the energy necessary for the meson to scatter from a nucleon at rest, with ran and tun the 
masses of the resonance and the nucleon respectively. Yr is the width of the resonance and gi is the isospin factor (2 
for isospin | and | for isospin | resonances). p(r) is the time dependent density of the system and appears due to 
application of the p-T-approximation, where T stands for the forward scattering amplitude. The important part for 
our purpose is the structure of the denominator, which represents a characteristic pole structure. We will replace to 2 
in the numerator by a constant factor since it causes divergence of Re£ rGt , only cured by application of subtracted 
dispersion relations corresponding to counterterms. The cutoff c also prevents divergences but leaving in the factor of 
u> 2 causes inaccurate numerical results. We show the spectral function for the case of the iV(1520) resonance, where 
fuNp ~ 7.0 |4fij | and the width of the resonance is T R = 120 GeV, in Fig. Elfor different densities at k = 0. 
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FIG. 2: Spectral function of the p-meson coupled to an N(1520)-hole pair for different densities. E = 582 MeVTn = 120 MeV 



We now discuss the system time r. Defining r = tl 1 2 t2 , as first done in [27|, seems to be a sensible choice, but 
when investigating the spectral function, one is dealing with retarded quantities and we do not want them to collect 
information from the future, i.e., we want to retain causality. This is why we choose r = t\ instead of the symmetric 
form t = tl t t2 , because in this case at a certain time r information of S ret that is located in the future of r enters 
the spectral function. This cannot happen with r = t\ as demonstrated in Appendix lAl 

For the self energy E < (r, ui) one should stick to the symmetric choice in order to fulfill the symmetry relation 
E < (ii, £2) = S > (t2, t\). However, doing this reduces accuracy in the numerical calculation, due to additional necessary 
Fourier transformations. Comparison of calculations using either T, < (ti, us) or £ < ( tl ^" t2 , u>) reveals minor differences 
in the rate when changing the temperature with system time t\ or tl ~^ t2 respectively, but on average, and hence in 
the final yield, the differences cancel, as could be expected. 



B. Contributions to the rate in time representation 



Before we calculate nonequilibrium rates and yields, we investigate how the present rate is created over time, i.e., 
we focus on which contributions to Eq. (|26|l come from which times in the past. Figs. 03 and 0] show the integrand 
for fixed energy uo = 750 MeV and fixed time r for relative times r — t. The first case shown (left of Fig. [3J is an 
equilibrium scenario with free p-mesons embedded in an environment at constant temperature. One can see that the 
main contribution comes from the very near past, but that there are also contributions from early times as well as large 
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oscillations leading to alternating positive and negative contributions. The second case (right of Fig. |3J) represents 
p-mesons, embedded in a constant heat bath, broadened and coupled to the N(1520) resonance at the time where the 
coupling and broadening was completely turned off again 10 fm/c ago in a way that will be discussed in Section ITlI CI 
(over a time of At = 7.18 fm/c, cf. Fig. 0- We use this duration because it results as the lifetime of the hadronic 
phase of the fireball from calculations described in Section UlI Fl and will hence be used in calculations of the dilepton 
yields). The integrand shows richer structures than in the equilibrium case, caused by the changing self energy in the 
past, which is quite remarkable: It shows that after having exposed the p-mesons to vacuum conditions for 10 fm/c, a 
memory of the situation in the more distant past is remaining. The third case (Fig. again represents free p-mesons 
but with the production in the heat bath, caused by , turned off 10 fm/c in the past. Comparison of the overall 
amplitude to that in Fig. [31 (left) shows that the rate at time r, resulting from the integral over the function shown 
in Fig. 0] (see Eq. will be very small compared to the rate before production was turned off, which is clear 

since the p-meson has a decay constant of 1.3 fm/c with that the rate is exponentially suppressed after production 
is turned off. The shape of the contribution is very interesting: Between the point where production was turned off 
and the present the integrand has a decaying (and oscillating) shape and contributions from the more distant past 
also remain. However, the largest values, contributing to the present rate, are situated close to the time, at which 
production was turned off. This is intuitively clear. The dileptons produced at the present time r most dominantly 
originate from decaying p-mesons, produced slightly before the time when production was turned off. Thereby the 
strong memory to that past time is visible in the integrand. 



~ 1-5 I 



> 
CD 

CD 



Q_ 

X 
CD 



1.0 



0.5 



0.0 



v -°- 5 
Q 

i -i.o 




m = 770 MeV 
r= 150 MeV 



co = 750 MeV 



past 



(present) 



10 

(x-t) [fm/c] 



15 



20 







> 

CD 

o 


0.3 - 
0.2 






H — ' 




1 

.3 


0.1 ■ 


exp( 


0.0 - 


-*— * 




V 

f- 
Q 


-0.1 - 


Re[i 


-0.2 - 




= 1520 MeV 
150 MeV 



= 240 MeV, r R2 =120 MeV 



m = 770 MeV, m, 
r, = 300 MeV, r, 

r R1 

co = 750 MeV 




past 



(present) 



10 

(x-t) [fm/c] 



15 



20 



FIG. 3: Contributions to the rate for a fixed frequency at given time r from times (t — t) in the past for the case of free p- mesons 
at constant temperature (left) and the case where the coupling to the N(1520) was turned off 10 fm/c before the present time 
(widths were also changed as indicated in the figure (see text) - index 1 refers to initial, 2 to final quantities) (right). 



A full interpretation of the various structures shown is difficult and it becomes clear that only the completely 
integrated yield represents a physical quantity. There is a dependence on the parameter A (cf. Eq. (|30() and related 
discussion) such that contributions from times further in the past are reduced using larger cutoffs A. On the other 
hand the accuracy of the numerics is strongly increased by employing larger cutoffs (when setting A to zero all the 
principal information is destroyed by diverging contributions, several orders of magnitude larger than the actual 
structure responsible for the creation of the rate). The early time contributions translate to diverging contributions 
for the smallest frequencies after Fourier transformation. The comparison of the rate from a free p-meson at constant 
temperature (i.e., the equilibrium case) calculated dynamically and the one calculated using the equilibrium formula 

dN . . 2e 4 m 4 1 
¥Z¥k (u) = (2^^" B( " MM ' (46) 

shows that the dynamically calculated rates, being integrals over strongly oscillating functions, as shown in Figs. 0and 
are reproduced very well by the numerics as nicely seen in Fig. [5] Here one also notices the large contributions at 
low frequencies, which cause the mentioned numerical uncertainties and therefore have to be reduced by introduction 
of the mentioned cutoff A. These contributions are being many orders of magnitude larger than those for frequencies 
at the upper end of the regarded range. 
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FIG. 4: Contributions to the rate for a fixed frequency at given time r from times (r — t) in the past. The production has been 
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FIG. 5: Comparison of the dynamically calculated rate (solid) to the one calculated from Eq. (1461 . the Markovian rate (dashed), 
shows very good agreement. 



C. Time scales of adaption for the spectral function, occupation number and dilepton rate 

In order to quantify the times that the mesons' spectral properties need to adjust to the evolving medium, we regard 
the cases of broadening and mass shifts and introduce a time dependent self energy that represents linear changes 
(see Fig. EJ in width or mass. We do not perform sudden changes using step functions in time because that would 
cause strong, hard to control oscillations in the Fourier transforms. Apart from the linear function, totally smooth 
functions like a hyperbolic tangent have been tested for the representation of the time evolution. It was found that 
the two kinks in the linear function do not cause additional problems when Fourier transforming. Hence, we chose 
this representation for the time evolution because it allows for the precise definition of the duration of the change, 
and by that of the time scale, which is not possible for absolutely smooth functions. 

The spectral function is given by the imaginary part of the Fourier transform in relative time coordinates of the 
retarded meson-propagator 

Df(ti,t 2 ) - Df(r,ti-t 2 ) -> -~ImDf(r,u) = A{t,u>), (47) 

where the system time r is chosen to be t\ for reasons discussed at the end of Section llH Al and in Appendix^] Taking 
t to be h-th leads to an about two times faster adjustment of the spectral function, which is caused by the fact that 
it is now able to gain information on the medium symmetrically from the past and the future, which is unphysical. 
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The spectral functions 's retardation with respect to medium changes is found by comparing the dynamically cal- 
culated spectral function at the time T g when the medium effects are fully turned off and the one calculated at the 
same time, assuming an instantaneous adaption. This spectral function is given directly by 



^Markov (t,U>) ImI}*M arkov 



1 



(r,u,) = - 



7T (oj 2 



Re5X Gt (r, u)f 



(48) 



The difference can be made explicit by calculating the difference in the moment J °° A(t b, ui) 2 L0 2 dLU of the two spectral 
functions or the difference in the peak position or height for mass shifts or broadening respectively. As an illustration of 
the procedure the spectral functions, which are compared, are shown for a particular example of in-medium broadening 
in Fig. [7| (full and dotted line). All methods lead to similar results. We find an exponentially decreasing difference 
with increasing duration of the change At (see Fig. EjJ. From this exponential drop we extract a time constant f, that 
for the case of broadening (and constant masses) is shown in Table[I]for different scenarios. Table[H]shows time scales 
for the adjustment of the spectral function to changes in the meson's mass, extracted from the peak position. Due to 
oscillations in the changing spectral function (see below) an exact extraction of a time scale is more complicated in 
this case, but the numbers still give a magnitude for the retardation. 

In the case of the p-meson (r vac « T2 — 150 MeV) we find a time scale of about 3 fm/c. That means that the 
behavior of the p mesons becomes adiabatic only for medium changes that are slow compared to the time of 3 fm/c, 
i.e., the spectral properties follow the changes in the medium nearly instantaneously only if the evolution is very 
slow as compared to the derived time scale. For real heavy ion reactions this means that the spectral properties of 
vector mesons will retain a certain memory of the past, and even if they decay outside the medium, they still carry 
a certain amount of information on the medium in that they were produced. This can become important, especially 
for w-mesons, having a vacuum width of 8.49 MeV [55|. Nevertheless, due to oscillations in the spectral function and 
rate and occurring interferences when calculating the dilepton yield (see below), this strong retardation may be put 
into perspective because the only physically meaningful influence of this retardation is that on the resulting yield (see 
We find the relation 



r cx — , with c between 2 and 3.5. 

12 



(49) 



c depends on the other quantities m and T± such that for a larger change in width, the spectral function takes 
longer to adjust to it. The reason for the time scale to be (at least) 2/T is that the spectral function is given by the 
imaginary part of the retarded propagator, which, representing a probability amplitude, is proportional to e~? Tt . Its 
square, an actual probability is proportional to e~ rt , giving the appropriate decay rate. 

This result can also be easily retrieved analytically when we assume a Breit-Wigner type self energy with constant 
width as well as a constant mass and solve the equation of motion that follows directly from Eq. I|34|l 



d tl D^(ti,t2)T 



(50) 
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FIG. 7: The two compared spectral functions right after the full change of the self energy is reached: dynamic (full) and 
Markovian (equals the final spectral function with F = 150 MeV) (dotted) calculation. The initial spectral function (dashed) 
is also shown. 
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TABLE I: Values of r for different masses, initial widths Ti and final widths IV 



by 

^(*i.*a) = -2©(*i - *a) 7= I e-^-^W^W-r^! - i a )) (51) 
H ' Vim 2 - T 2 2 

and we find the decay width to be T/2. The numerical solution shows the same behavior, but here we can go further 
and apply changes e.g. to the mass and see how the propagator behaves in time. In Fig. [8] we show the numerical 
solution for the retarded propagator in time representation D™^(r, t) at different fixed times r, for a constant width 
r = 150 MeV and a mass, changed from 400 to 770 MeV (n before the change of duration At = 7.18 fm/c started, 
T2 shortly (1 fm/c) after the finished change, and T3 long (ss 75 fm/c) after the change). We combine the plots of this 
propagator with an exponential -^=^=p e~? rt to show that it decays with the analytically found decay constant. 
One can also nicely see how the wavelength, related to the mass by A = ~^|=f r 2 ( see Q^lflh changes with time 
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TABLE II: Values of r for different widths and mass shifts Am to masses below the final vacuum mass 7712 = 770 MeV. 
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FIG. 8: Numerical results for the retarded meson propagator in the two-time representation for a constant mass of 400 MeV, 
shortly (1 fm/c) after a mass shift to 770 MeV, and a rather long time (about 75 fm/c) after that mass shift. Discussion in the 
text. 



and that at earlier times in the second plot (j2 — t large) the frequency still corresponds to masses close to the initial 
mass and at times close to the present to the final mass. This is a nice illustration of the memory of the retarded 
propagator and hence the spectral function in the two-time representation. The stronger retardation (c > 2) for larger 
changes in the self energy can be explained by considering that the change of the width or mass introduces additional 
time scales that might as well become of importance such that it is no longer so easy to extract the time scale directly 
and only from the width. 

For the occupation number of the vector mesons, which is basically given by its propagator G < (t, lu) and has to 
be compared to Gj^ arkov (t, u>) — 2inB{T(r), w ) ^.Markov (V, lu) at time r fj, we find a by about 50% faster adaption to 
medium changes: The quantity is not retarded and gains information from the future when Fourier transforming. On 
the other hand, for the causal dilepton rate (|26|) . which is compared to Eq. I|46(l with the Bose distribution ub(T(t), lu) 
and the spectral function ^Markov (t, lu) at time r ff, this is not the case: We find a retardation much alike that of the 
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spectral function with similar time scales as given in Eq. (|49ll . 



D. Quantum interference 



As already mentioned, there occurs another interesting effect, arising from the quantum mechanical nature of 
nonequilibrium dilepton emission. Oscillations in the changing spectral functions, occupation numbers and production 
rates appear as well as interferences that cannot be present in an approximate semi-classical calculation. The fact 
that negative values occur temporarily is comparable to the well known observation that the Wigner function, the 
quantum mechanical analogue to the classical phase space distribution, is not necessarily positive definite (see e.g. 
|56j). It is interesting to note that although the rates may oscillate below zero, the total yield will stay positive, being 
the only physical observable calculated. It can be shown that the accumulated yield is indeed proportional to the 
square of an amplitude [30l | . The rate in this case is not to be confused with the rate that results from measuring the 
yield and simply dividing by the time interval. Instead, the rate (|26ll calculated here has the full quantum mechanical 
information incorporated and contains interferences that can cause cancellations - in fact, the rate has to be able to 
become negative. An example for the occurring oscillations is depicted in Fig. [!5]for the spectral function and the rate 
for the case of a mass shift to 400 MeV in the medium at a constant temperature of 175 MeV. The resulting yield 
integrated over the time interval beginning with the onset of the changes and ending at the indicated time is shown 
in Fig. ED It can be nicely seen that the yield stays positive throughout the evolution while the temporary rate can 
become negative. 
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FIG. 9: Spectral function and rate for the change of the mass from m = 400 MeV to 770 MeV (constant r=150 MeV and 
constant T = 175 MeV) directly after the self energy has reached its final form (after 7.18 fm/c) and 2 fm/c later. Oscillations 
and negative values appear in the intermediate spectral functions and rates. 



E. Yields at constant temperature 



In order to get a first impression of how memory effects can affect dilepton yields, we investigate yields from 
systems with a constant size at constant temperature. This helps to understand the more complex situation of an 
evolving fireball, treated in the next section. For five different scenarios we perform comparisons of the dynamic to 
the Markovian calculations in that the spectral properties adjust instantaneously to the medium and the dilepton 
rate has no memory: Modification of the p-width to 400 MeV in the medium, shift of the p-mass using a constant 
coupling at the p-7*-vertex to 400 MeV in-medium mass, coupling of the p to the N(1520) resonance with and without 
additional broadening of the resonance and the p-meson. For the w-meson we consider a mass-shift by 100 MeV to 
682 MeV and broadening to 40 MeV 52]. We note that the situation for the 0- meson might be similar. In each case 
the vacuum spectral function is approached linearly within the interval of 7.18 fm/c, over that we integrate the rates. 
The temperature is set to 175 MeV. The results for the first four (p-meson-) scenarios are shown in Fig. ^2 The 
largest, most noticable difference is found for the mass shift of the p-meson: The yield from the dynamic calculation 
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FIG. 10: Yields corresponding to rates from Fig. [5] (and some more) integrated over time intervals starting with beginning 
changes and ending after the indicated duration (3.18 fm/c, 5.18 fm/c, 7.18 fm/c, 9.18 fm/c, 11.18 fm/c and 15.18 fm/c). 



is increased by a factor of about 1.8 in the range from 200 to 450 MeV due to the inherent memory integrals. This 
can be understood since as opposed to the Markov case the spectral function has a certain memory of the in-medium 
conditions and the stronger weighting of lower masses due to the Bose factor increases the effect additionally. 

The other cases show differences, but not as pronounced. We see a minor enhancement of the lower mass tail 
due to the Bose factor in all cases. This is because for the dynamic case there remains a memory to the larger 
in-medium width and the resonance peak, giving more strength to the more strongly weighted lower mass part of 
the distribution than in the Markovian case. In the case of strongly broadened p-mesons, one can nicely see that the 
dynamically calculated yield possesses a broader distribution and that for the case of the coupling to the N(1520) 
without broadening the resonance peak is stronger (by about a factor of 1.5), due to the system's memory to the 
coupling in the medium. Here, in the dynamic calculation the vacuum peak is also shifted further to the right due to 
memory of the level repulsion effect. For the case of coupling to the N(1520) with additional broadening this effect 
is harder to see because of the very broad distribution - however, the yield around the resonance peak is stronger in 
this case as well. 

For the w-meson (Fig. 1121) large differences show up due to the strong retardation related to the small widths 
involved. One nicely sees that the in-medium peak is more pronounced in the yield from the dynamic calculation with 
memory, whereas the final peak is not yet visible. The stronger yield above 800 MeV can be explained by the memory 
to the larger in-medium width, whereas the decreased yield between 450 and 650 MeV must be due to interferences, 
discussed in Section IlII Dl The Bose factor does not have that much influence in this case because the rather small 
mass shift and broadening do not produce a spectral function with much weight at small frequencies. 

To conclude this section we state that for the p-meson we indeed found noticeable modifications of the yields that 
one would expect qualitatively when including finite memory for the spectral properties. For the w-nieson, the much 
longer time scales for the adjustment to the medium further cause structures resulting from quantum mechanical 
interference after a comparably quick change over 7.18 fm/c. 



After having had a first exploration of the influence of finite memory on the dilepton production for constant 
temperature, we now study a more realistic situation in that memory effects are expected to be of importance. The 
hadronic phase of a fireball, created in a heavy ion collision, lives for about Tfi ro b a ii = 5 to 10 fm/c at SpS energies. 
This is the time in that the changes of the vector mesons' self energy take place, a time comparable to the derived 
time scale of retardation for the spectral properties and hence the dilepton rate (see Section lill CI) : 



F. Fireball model and dilepton yields 



Tfiroball — T 



(52) 



This is why the consideration of memory effects is important for heavy ion collisions. 
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FIG. 11: Comparison of the dynamically calculated (solid) to the Markovian (dashed) dilepton yields from an interval of 
duration At = 7.18 fm/c in that the self energy was changed linearly from initial to final values given in the figure. Index 1 
indicates initial values, index 2 final ones. The temperature was kept constant at T — 175 MeV. 



We will model the fireball evolution and convolute it with the calculated time dependent rates. We choose for the 
effective volume a longitudinal Bjorken expansion combined with an accelerating radial flow 

Kfffr > t ) = ircT(r Q + v q (t - t ) + 0.5a (T - t ) 2 ) 2 , (53) 



with To = 6.5 fm, vo = 0.15 c and ao = 0.05 c 2 /fm (see also |57ll58) ). From Q53I1 and the constraint of conserved entropy 
(given by a constant entropy per baryon S/A = 26 for SPS energies [5^. l6dL l6l| ) . temperature and chemical potentials 
follow as functions of time. The initial temperature of the fireball is taken slightly above that at chemical freezeout, 
175 MeV, whereas the final temperature, reached after a lifetime of about 7.18 fm/c is 120 MeV (thermal freezeout). 
At this point, we turn off further dilepton production, by decreasing (quasi instantaneously) the temperature towards 
zero. Afterwards only existing vector mesons decay. We take Ab — 70 as the number of participant baryons per 
unit rapidity and include the 56 lightest baryonic and mesonic states in this calculation. Using the calculated time 
dependent temperature T(t) in the computation of the rate and the given volume T4ff(T), we can integrate the rate 
and get the accumulated yield per unit four momentum. 

Before this we want to consider the dynamic process of decay, occuring after the production is turned off. Therefore 
we prepare a situation of constant temperature and fixed spectral properties for the p-meson and turn down the 
temperature at a certain point. We then compare the resulting final yield to that calculated analytically for the 
Markov case using 

dN , . dN , , / 1 \ 2e 4 mi 1 



The result is given in Fig. ^| It nicely shows the difference in the decay process: The dynamic calculation leads to 
a much more pronounced peak while the Markov approximation has larger contributions away from the peak. This 
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FIG. 12: Comparison of the dynamically calculated (solid) to the Markovian (dashed) dilepton yield from oj-mesons from an 
interval of duration At — 7.18 fm/c in that the self energy was changed linearly from initial to final values given in the figure. 
Index 1 indicates initial values, index 2 final ones. The temperature was kept constant at T = 175 MeV. 



can be understood by the investigating the dynamic evolution of the occupation of p- mesons, given by D<(t,u>), after 
turning down the temperature. We show this in Fig. 1141 It nicely exhibits how the "thermal" occupation far from 
the pole mass disappears. This is not given in the approximate Markov formula (Eq. 1)54(1 ). where ?t-b(w) is constant 
in time. Integrating the yield over invariant masses M from 125 to 1500 MeV yields the same overall dilepton number 
for both cases within a percent. Furthermore we find the expected behavior of the dilepton yield integrated over the 
invariant mass after the time of switching off r ff: Its increase per unit time decreases with e~ Tt (for t > r tf). 
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FIG. 13: For an equilibrium situation (p-meson vacuum 
spectral function and constant temperature) after turning 
down the temperature, the dynamically calculated yield 
(solid line) differs from the analytically calculated one 
(Eq. 1541 ') (dashed line). Also shown is the yield at time 
r ff from an interval starting 7.18 fm/c before (same for 
both calculations) (dotted line). 



FIG. 14: Evolution of D^-(t,u) after having turned off 
the production at time T g. 



Now we turn to the situation of an expanding fireball. We consider five different scenarios as in the constant 
temperature and volume case (see Section fill E|) : Broadening of the p- meson to a width of 400 MeV, a mass-shift 
of the p-meson to 400 MeV mass, coupling of the p-meson to the N(1520) resonance with and without broadening, 
and modification of mass and width of the w-meson. Comparison to calculations at constant temperature reveals the 
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additional influence of dropping temperature and increasing volume, as well as that of the mentioned difference in 
behavior to the Markov calculation after the production has been turned off. Overall the differences in the yields are 
stronger for the case of changing temperature and volume. Fig. El shows the comparison to the Markov calculation 
for the four p-meson scenarios. There always remains a certain memory of the higher temperatures at earlier times 
which together with the increasing volume leads to an increased yield within the intermediate mass regime from 400 
to 700 MeV for the case of a mass shift to 400 MeV in-medium mass, which was mostly restricted to the vicinity of 400 
MeV in the constant temperature case. The dynamically calculated yield around the in-medium peak is now a factor 
of two larger than that from the Markovian calculation. This constitutes a significant effect and shows the importance 
of the consideration of the full dynamics for mass shifts, especially if one asks for precise theoretical predictions. In 
the case of coupling to the N(1520) resonance with and without broadening one now gets more pronounced s-shape 
structures in the yield than in the Markov case. In addition the different behavior after the temperature was turned 
off becomes visible in the peak at the vacuum mass. This effect even dominates for the case of the broadened p-meson 
such that the actual effect of a broader distribution in the dynamic case, which was seen in the constant temperature 
calculation, is hidden. 

For the situation of a propagating w-meson mode (Fig. we calculate the final Markovian yield by adding the 
contribution given in Eq. I|54|) but using the initial temperature of 175 MeV, since most of the states, populated 
at high temperature, are still decaying at the time of thermal freeze-out and contribute to the final yield. One can 
see an enhancement by a factor of more than two between the in-medium mass and about 750 MeV in the dynamic 
calculation, caused by the finite memory to the in-medium properties. Due to the finite resolution, the numerical 
treatment can only approximate the narrow u vacuum peak and gives visible wiggles in the yield. 




FIG. 15: Comparison of the dynamically calculated (solid) to the Markovian (dashed) final dilepton yields where the self energy 
was changed linearly over an interval of duration At — 7.18 fm/c as indicated in the corresponding figure (index 1 for initial, 
2 for final quantities) . The time dependent rate was then convoluted with the temperature and volume of the fireball. Dotted 
lines show yields from unmodified (vacuum) p-mesons. 



What about the non-zero momentum modes? The shown calculations for the momentum mode k — can be 
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FIG. 16: Comparison of the dynamically calculated (solid) to the Markovian (dashed) final dilepton yields for the tj-meson. 
The self energy was changed linearly over an interval of duration At = 7.18 fm/c as indicated in the figure (index 1 for initial, 
2 for final quantities). The time dependent rate was then convoluted with the temperature and volume of the fireball. 



extrapolated at least to further lower momentum modes, since the self energy is a continuous function with respect 
to momentum. Typically, for the lowest momenta the strongest medium modifications are seen in experiment [3, 62] . 
The yield per unit rapidity and invariant mass can then be given by converting variables and approximating the 
integral over transverse momenta by a product of the yield for momentum k = with the regarded momentum range 
Ak_L. In order to get the yield for y = in the center of mass frame, in the Bjorken model one has to integrate 
over all rapidities of the fluid cells, using that y = r) + y', where rj is the space-time rapidity of the cell and y' the 
momentum rapidity of the particles in the restframe of the cell. 



dN 



dMdy 



2; 



Ml dn I dy' I d 2 k ± (-^-) 8{r, + y') 

k±<k cut ,y=0,CM J-oo J-oo JO \ a£,a K J LR.fey {y'),k ± 

r°° r°° f dN \ 
«M/ d v / dy' — 3- S( V + y')irk 2 cnt (55) 

J-oo J-oo \ a£jCL K / LR : k n (y'),k ± =0 



where CM indicates the center of mass frame, LR the local restframe and fc cu t is the upper momentum cutoff for 
transverse momenta. The cutoff k cut is chosen to be 500 MeV in 0, for investigation of dilepton yields for lower 
transverse momenta. For the Bjorken model we encounter the problem that arbitrarily large rapidities y' and hence 
also large longitudinal momenta contribute. Approximating the yields for these momentum modes by that for k = 
is not reasonable. 

On the other hand for a "semi-static" expanding fireball model as employed in e.g. 0,0] the situation becomes 
much simpler as in this case the only dilepton longitudinal momentum contributing to the midrapidity yield is fc|| = 
and we have 



dN 



dM dy 



„ f dN \ 

M / d 2 k ± ^ 

Jo ^{dEd^k), 



fcj_<fccut,J/=0 J vwaju, n, / k = Q,fc 



M (^A (56) 



\dEd 3 k J 



k„=0,k i =0 



No momenta larger than fc cut are involved for this case and one can approximate the yields for the appearing small 
momenta by those for zero momentum. Hence also for measurable finite momentum intervals the resulting modifi- 
cations of the yields due to finite memory are expected to be completely similar and our primary message remains 
valid. 

We summarize by stating that inclusion of the fireball evolution leads to further enhancement of the differences 
between the dilepton yields calculated using a full noncquilibrium quantum field theoretical formulation including 
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retardation and its much simpler Markovian approximation. For the realistic situation of a heavy ion collision memory 
effects have to be considered when it comes to predicting medium modifications of particles. 

IV. SUMMARY AND CONCLUSIONS 

In the present work, we have numerically calculated dilepton production rates within a nonequilibrium field theory 
formalism, based on the real time approach of Schwinger and Keldysh |23j, 124. \2aL I2(t| . We employed the Kadanoff- 
Baym equations, generalized to the relativistic Dirac structure, to derive the formula for the dilepton production 
rate, which is non-local in time and hence includes the usually disregarded finite memory. The rate is in principle 
a (half) Fourier transform over past times of the virtual photon occupation, described by the Green function D< in 
the two-time representation. Medium modifications of the vector mesons enter this production rate via the virtual 
photon's self energy that is connected to the vector meson propagator by the principle of vector meson dominance. 

The equation of motion for the retarded and advanced vector meson propagator also followed from the Kadanoff- 
Baym equations. From its solution it was possible to extract the dynamic behavior of the meson's spectral function 
and to define a time scale on that it adjusts to medium changes, which were incorporated by the definition of a time 
dependent self energy for the vector mesons. This time scale was found to be proportional to the inverse vacuum 
width of the meson like c/F. c lies between 2 and about 3.5. The time that the dilepton rate needs to follow 
changes was found to be approximately equal to that for the spectral function. Since these time scales lie in the 
range of typical times in that the hadronic phase of fireballs in heavy ion reactions exists and evolves, we expected an 
influence of the found retardation on the dilepton yields and investigated different scenarios quantitatively. Possible 
medium modifications of the vector mesons, as shifted pole masses, motivated by Brown-Rho scaling, broadening and 
coupling to resonance-hole pairs were considered. Investigation of the dynamical off-shell evolution of the spectral 
function, meson occupation number and resulting dilepton rate revealed the quantum mechanical nature of the system, 
manifested in appearing oscillations and hence interferences in the regarded quantities after changes of the self energy 
that were fast as compared to the derived time scales. All quantities were found to possibly become negative, which 
is a major difference to the Markov approximation within which all quantities are always positive definite. The 
oscillations potentially cancel when the rate is integrated over time such that the measurable dilepton yield is found 
to be positive, as it has to be. On the other hand, for situations where the particular rate can become negative, a 
semi-classical interpretation with positive definite rates, as employed in present day transport codes, is not possible. 

We first calculated dilepton yields from a time dependent system at constant volume and temperature considering 
different possible medium modifications and compared to the yields found assuming an instantaneous (Markovian) 
adjustment of all quantities to the medium. For the p-meson we found qualitatively expected results for each case. 
Quantitatively, the strongest difference in the differently calculated yields appeared for the case of a mass shift to 400 
MeV in-medium mass, as suggested by Brown-Rho scaling. This is due to the stronger weighting of the lower masses 
by the Bose factor and since the spectral function moves more slowly towards the free mass in the dynamic situation, 
we found more enhancement in the yield than for the Markov calculation without memory. It was shown that there 
can be more than a factor of two difference. For the medium modified w-meson the qualitatively expected results plus 
additional structures caused by quantum interferences were seen. These are stronger than in the case of the p-meson 
because the performed change is faster relative to the time scale for adiabatic behavior of the u-meson, which has a 
very small width. 

The introduction of a fireball model allowed us to perform more realistic simulations including changing temperature 
and volume, adjusted to the case of the SpS with energies of 158 AGeV. The modifications of the yields due to inclusion 
of finite memory were enhanced as compared to the constant temperature case. Also the behavior of the regarded 
quantities after freeze-out was revealed to be different from that assumed in Markovian calculations. For the mass 
shift of the p-meson a factor of two difference in the two calculated yields was found around the in-medium mass - a 
significant effect. For the case of substantial broadening without further modifications, the differences are only very 
small. In all scenarios, also for the w-meson, lower invariant masses were slightly enhanced in the dynamic calculation. 
For the coupling to the N(1520) resonance the peaks were more pronounced than in the Markov calculation. 

In summary, our findings show that exact treatment of medium modifications of particles in relativistic heavy ion 
collisions requires the consideration of memory effects. This is in particular true for mass shifts or the occurrence of 
two-peak structures in their spectral properties. It has to be a future task to find how to incorporate such memory 
effects in semi-classical transport simulations. This should be also of relevance for the description of production and 
propagation of vector mesons through cold nuclei in photo-nuclear reactions [63l Ir34| . 
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APPENDIX A: ON THE ACAUSALITY OF THE KADANOFF-BAYM PARAMETRIZATION 

In this appendix we first show that the choice of time variables tkb = "t ^ and At = t a — tp leads to an acausal 
spectral function A(tkb,(j). This choice will be referred to as the Kadanoff-Baym parametrization since it was first 
introduced by Kadanoff and Baym in 27] for a first order gradient expansion of the full quantum transport equations. 
We introduce a different parametrization r = t a and At = t a ~ tp and obtain a spectral function A(t, u) that has no 
information on the future incorporated. 

The following analysis of the differential equation of the retarded propagator 



D tA D?(tA,t B ) 



dt 



2 Zf(t A ,t 2 )Df(t 2 ,t B ) 



(Al) 



allows us to find the times, from which medium information T, ret (ti,t 2 ) is contributed to D p ct (tA, t B ). 



D r p et (tA,t B ) follows by integration of D tl D T p ct (ti,t B ) = di 2 £p et (*i; ^i^'fe, is) over tie[t B ,t A ], starting at 
t\ = t B 2 . Since all quantities are retarded, we always have t A > t\ > t 2 > t B . Z? p ct (ii,is) itself has contributions 
from -Dp Ct (t2,^s) and Y7 p ct (ti,t 2 ) with t 2 e[t B ,ti]. Fig. El schematically shows where in time these contributions are 
located. On the horizontal line enter contributions from £ p et (ii, t 2 ), on the vertical line those from Z) p ct (i2, t B ). After 
integration over tie[t B ,t A ], one finds that D^^a^s) is created solely with information T, p ot (ti,t 2 ), defined at times 
lying in the triangle shown in Fig. 1181 



t A 
t B 




Df(t A ,t B ) + - 
Df(t u t B )± y 

Df{h,t B ) k 

i / D ret 

p 




FIG. 17: (Color online) Ranges in time (shaded regions) 
that contribute to D p ot (t\,ts) in the t^-i^-plane, where 
ta {tp) stands for the first (second) time argument. 



FIG. 18: (Color online) Contributions of E rot (U, tg) to 
Dp et (tA,tB) come from times t a , tp, lying within the 
shaded region. 



When we now choose to transform time variables like tkb = t ° L ^ 2 tf3 and At = t a — tp, we find that the Green 
function D p et (rKB, At) is created also with information defined in its future as demonstrated in Fig. El Hence, 
this parametrization cannot be causal and is not suitable for finding appropriate time scales for adaption of the 
spectral function to medium changes in time. The spectral function A{t-ks,uj) = ~^I m -Dp ct ( r KB, ^) is the imaginary 
part of the Fourier transform in relative time coordinates of the Green function and hence would have incorporated 
information on self energies from the future. 



2 For the case of Dt 1 = (-9?, — m2 ) the initial conditions are D p et (tg,tg) = and dt 1 -D p et (ti, tg)\t 1 ^t B = — 1- 
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FIG. 19: (Color online) After transformation to tkb = 
f °+^ anc j At — t a — tp the Green function at time tkb = 
tA+ 2 tB , D 1 ^ (r KB , At), has contributions of E rot (t Q , tp) 
also from its future, indicated by the upper shaded region. 



FIG. 20: (Color online) Choosing the parametrization 
t — t a and At = t a — tp leads to a Green function 
Dp ct (r = t A ,At) without contributions of E rct (t a , tp) 
from the future. 



On the other hand, using r = t a and At = t a — tp, as done throughout this work, the situation looks as depicted 
in Fig. 1201 Only past times contribute to the Green function at time t — tA- A spectral function defined as 
A(t a ,tj) = — ilm_Dp Ct (t Q , lu) is causal as required. 
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